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5 ■ Abstract 
CN 

£> \ We describe how Grobner bases can be used to solve the reduction 

problem for Feynman integrals, i.e. to construct an algorithm that 
provides the possibility to express a Feynman integral of a given family 
as a linear combination of some master integrals. Our approach is based 
ly-j \ on a generalized Buchberger algorithm for constructing Grobner-type 

O \ bases associated with polynomials of shift operators. We illustrate it 

through various examples of reduction problems for families of one- and 
two-loop Feynman integrals. We also solve the reduction problem for a 
^h, family of integrals contributing to the three-loop static quark potential. 
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1 Introduction 



The important mathematical problem of evaluating Feynman integrals arises natu- 
rally in elementary-particle physics when one treats quantum-theoretical amplitudes 
in the framework of perturbation theory. This problem originated in the early days 
of perturbative quantum field theory. Over more than five decades, a great variety 
of methods for evaluating Feynman integrals has been developed. However, to check 
whether the Standard Model or its extensions describe adequately particle interac- 
tions observed in experiments, one needs to perform more and more sophisticated 
calculations, so that one tries not only to update existing methods but also develop 
new effective methods of evaluating Feynman integrals. 

After a tensor reduction based on some projectors a given Feynman graph gener- 
ates various scalar Feynman integrals that have the same structure of the integrand 
with various distributions of powers of propagators which we shall also call indices. 
Let F(ai, a,2, ■ ■ ■ , a n ) be a scalar dimensionally regularized pQ Feynman integral cor- 
responding to a given graph and labelled by the (integer) indices, af. 



where ki, i = 1, . . . , h, are loop momenta and the denominators are given by 



with r = l,...,n. The matrix A l J depends on the choice of the loop momenta. 
The momenta pi are either the loop momenta pi — ki, i = 1, . . . , h, or independent 
external momenta Ph+i — Qi, ■ ■ ■ ,Ph+n = q n of the graph. Irreducible polynomials in 
the numerator can be represented as denominators raised to negative powers. For 
example, the denominator corresponding to the propagator of a massless particle is 
k 2 = k^ — k 2 . Usual prescriptions k 2 = k 2 +i0, etc. are implied. Formally, dimensional 
regularization pQ is denoted by the change d 4 k = dk^dk — > d d k, where d = 4 — 2e is 
a general complex number. The Feynman integrals are functions of the masses, mi, 
and kinematic invariants, ■ qj. However, we shall omit this dependence because we 
shall pay special attention to the dependence on the indices. We shall also omit the 
dependence on d. 

A straightforward strategy is to evaluate, by some methods, every scalar Feynman 
integral resulting from the given graph. If the number of these integrals is small this 
strategy is quite reasonable. In non-trivial situations, where the number of different 
scalar integrals can be at the level of hundreds and thousands, this strategy looks too 
complicated. A well-known optimal strategy here is to derive, without calculation, 
and then apply some relations between the given family of Feynman integrals as 
recurrence relations. A well-known standard way to obtain such relations is provided 




(1) 




(2) 



i>j>l 



1 



by the method of integration by parts (IBP) [2] which is based on the fact that any 
dimensionally regularized integral of the form 



is equal to zero. Here / is the integrand in (JI|) . More precisely, one tries to use 
IBP relations in order to express a general dimensionally regularized integral of the 
given family as a linear combination of some 'irreducible' integrals which are also 
called master integrals. Therefore the whole problem decomposes into two parts: 
the construction of a reduction algorithm and the evaluation of the master Feynman 
integrals. 

There were several recent attempts to make the reduction procedure systematic: 

(i) Using the fact that the total number of IBP equations grows faster than the 
number of independent Feynman integrals, when one increases the total power of the 
numerator and denominator, one can sooner or later obtain an overdetermined system 
of equations E] which can be solved. (There is a public version of implementing 
the corresponding algorithm on a computer 

(ii) Using relations that can be obtained by tricks with shifting dimension jB]. 
(Hi) Baikov's method [7j. 

Another attempt in this direction is based on the use of Grobner bases j^]. The 
first attempt to apply the theory of Grobner bases in the reduction problems for 
Feynman integrals was made in ^U|, where IBP relations were reduced to differential 
equations. To do this, it is assumed that there is a non-zero mass for each line. For 
differential equations one can then apply some standard algorithms for constructing 
corresponding Grobner bases. 

In [TIH ITT] it was pointed out that the straightforward implementation of the 
Buchberger algorithm in the case of IBP relations is problematic because it requires 
a lot of computer time even in simple examples. One of the possible modifications 
is related to the Janet bases |T2]. We are going to modify the Buchberger algorithm 
in another way, taking into account explicitly such properties as boundary conditions 
(which characterize all the regions of indices where the Feynman integrals are equal 
to zero), so that it will be possible to apply it to solve the reduction problem in 
complicated situations. 

In the next section, we shall briefly describe what the Grobner basis and the Buch- 
berger algorithm are in the classical problem related to solving systems of algebraic 
equations. In Section 3, we shall turn to IBP relations and describe our strategy of 
constructing Grobner bases associated with the given problem, with the help of a 
modification of the standard Buchberger algorithm. We shall explain how these re- 
sults can be applied to solve IBP relations. We shall illustrate our strategy, through 
various examples, in Section 4. In particular, we shall apply our algorithm to a family 
of three-loop Feynman integrals with a one-loop insertion relevant to the three-loop 
quark potential. We shall also evaluate the master integrals using the method based 
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on Mellin-Barnes representation. In Conclusion, we shall characterize the status and 
perspectives of our method. 

2 Grobner basis and Buchberger algorithm 

The notion of the Grobner basis was invented by Buchberger [S] when he constructed 
an algorithm to answer certain questions on the structure of ideals of polynomial 
rings. 

Let A = C[xi, . . . ,x n ] be the commutative ring of polynomials of n variables 
xi, . . . ,x n over C and X C A be an ideal 3 . A classical problem is to construct an 
algorithm that shows whether a given element g G A is a member of X or not. A 
finite set of polynomials in X is said to be a basis of X if any element of X can be 
represented as a linear combination of its elements, where the coefficients are some 
elements of A. Let us fix a basis {fx, f 2 , ■ ■ ■ , fk} of X. The problem is to find out 
whether there are polynomials r\, . . . , G A such that g = rxfx + . . . + r^f^. 

Let n = 1. In this case any ideal is generated by one element / = a + a\x + 
a 2 x 2 + . . . + a m x m . Now if we want to find out whether an element g = b + bxx + 
b 2 x 2 + . . . + b\x l can be represented as rf we first check if I > m. If so, we replace 
g with g — (bi/a m )x l ~ m f , 'killing' the leading term of g. This procedure is repeated 
until the degree of a 'current' polynomial obtained from g becomes less than m. It is 
clear that the resulting polynomial is equal to zero if and only if g can be represented 
as rf. 

Now let n > 1. Let us consider an algorithm that will answer this problem for 
some bases of the ideal. (We will see later that this problem can be solved if we 
have a so-called Grobner basis at hand.) To describe it, one needs the notion of an 
ordering of monomials cx 1 ^ . . . x % ™ where c G C and the notion of the leading term 
(an analogue of the intuitive one in the case n — 1). In the simplest variant of 
lexicographical ordering, a set (ii, . . . , i n ) is said to be higher than a set (ji, . . . , j n ) if 
there is I < n such that i\ = ji, i 2 — ji-, ■ ■ ■ , k-i = jl-i an d h > ji- The ordering 
is denoted as (ii, . . . , i n ) y (ji, ■ ■ ■ ,j n ). We shall also say that the corresponding 
monomial cx 1 ^ . . . x l ™ is higher than the monomial dx 3 ^ . . . x^ . 

One can introduce various orderings, for example, the degree-lexicographical 
ordering, where (zi,...,z n ) >- (j 1 ,...,j n ) if > 12 3k, or £z fc = and 

(ii, . . . ,i n ) y (ji,...,j n ) in the sense of the lexicographical ordering. The only two 
axioms that the ordering has to satisfy are that 1 is the only minimal element under 
this ordering and that if fi y f 2 then gf t y gf 2 for any g. 

Let us fix an ordering. The leading term (under this ordering) of a polynomial 

P{x\, . . . , Xn) = ^ Ci lr .. t i n Xi . . . X n 

3 A non-empty subset I of a ring R is called a left (right) ideal if (?) for any a, b 6 X one has 
a + b G X and (ii) for any a G X,c G R one has ca G X (ac G X respectively). In the case of 
commutative rings there is no difference between left and right ideals. 



3 



is the non-zero monomial such that the degree . . . , i°) is higher 

than the degrees of other monomials in P. Let us denote it by P. We have P = P+P, 
where P is the sum of the remaining terms. 

Let us return to the problem formulated above. Suppose that the leading term 
of the given polynomial g is divisible by the leading term or some polynomial of the 
basis, i.e. g = Qfi where Q is a monomial. Let g\ = g — Qfi- It is clear that the 
leading term of gi is lower than the leading term of g and that g\ G X if and only if 
g G X. One can go further and proceed with g\ as with g, using the same f\ or some 
other element fj of the initial basis, and obtain similarly g 2 ,g^,.... The procedure is 
repeated until one obtains g% = or an element g\ such that g\ is not divisible by any 
leading term fa. We will say that g is reduced to gi modulo the basis f 2 , ■ ■ ■ , fk}. 

A basis {fa, f 2 , . . . , fj.} is called a Grobner basis of the given ideal if any polynomial 
g G X is reduced by the described procedure to zero for any sequence of reductions. 
Given a Grobner basis we obtain an algorithm to verify whether an element g G A 
is a member of X. There are many other questions on the structure of the ideal that 
can be answered constructively if one has a Grobner basis, but they are beyond the 
topic of the paper. 

Generally a basis is not a Grobner basis. Let fa = x\ and fa — 1 + x\ and let 
X be generated by fa and f 2 . It is easy to verify that {fa, fa} is a Grobner basis of 
X. Now let f[ = x\x 2 . The set {f[, f 2 } is again a basis of X (indeed, f[ = x 2 fi and 
fi = —x 2 f[ + x\f 2 ). However, {f[,f 2 } is not a Grobner basis because the element 
X\ G X cannot be reduced modulo {/{, / 2 }. 

On the other hand, given any initial basis {fi, f 2 , ■ ■ ■ , fk} of the ideal X one can 
construct a Grobner basis starting from it and using the Buchberger algorithm which 
consists of the following steps. 

Suppose that /, = wqi and fj = wqj where w, qi and qj are monomials and w is 
not a constant. Define /y = /jg^ — /^gj. Reduce this polynomial modulo the set {fi} 
as described above. If one obtains a non-zero polynomial by this reduction, add it to 
the set {fi}. Consider then the other elements with f[ = wq[ and /j = w'q'j for some 
non-constant w' . If there is nothing to do according to this procedure one obtains 
a Grobner basis. It has been proven by Buchberger that such a procedure stops 
after a finite number of steps. 

The Buchberger algorithm can take much computer time to construct a Grobner 
basis, but once it has been constructed, one can use the reduction procedure which 
works generally much faster. 

To conclude, the problem formulated in the beginning of this section can be solved 
by choosing an ordering and constructing the corresponding Grobner basis using the 
Buchberger algorithm. After that, one applies the reduction procedure modulo the 
constructed Grobner basis to verify whether a given element belongs to the given 
ideal X. 
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3 Reduction problem for Feynman integrals 



Practically, one uses relations (JHJ) of the following form: 

d ^...^( B _l_) =0 . (4) 

Here E r are denominators in (fl)l. ki, . . . ,kh are loop momenta and p\ — k±, . . . ,Ph — 
kh,Ph+i — Qi, ■ ■ ■ ,Ph+n = In, where q±, . . . ,q n are independent external momenta. 

After the differentiation, resulting scalar products, ki • kj and ki ■ qj are expressed 
in terms of the factors in the denominator, by inverting (j2J), and one arrives at IBP 
relations which can be written as 

^2 Ci F {o-i + ■ ■ ■ , «n + &i, n ) = , (5) 

where 6, i are integ polynomials in a,-, d, masses and kinematic invariants, 

and F(a\, . . . , a n ) are Feynman integrals of the given family. These relations can 
be written in terms of shift operators i + and i~ which are defined as 

i ■ F(a 1 , a 2 , . . . , a n ) = F(a 1 , . . . , a^i, ai ± 1, a i+1 , . . . , a n ) . 

At this point, we would like to turn from the 'physical' shift operators to 
'mathematical' shift operators. (We believe that the physical notation can be am- 
biguous: for example, it is not immediately clear whether the operators are applied 
to a function of the indices, or to some of its values.) 

Let /C be the field of rational functions of physical variables m^, q^ ■ qj, d, and 
A be the algebra 4 over K, generated by elements Yj, Y~ and Ai with the following 
relations: 

^i^j = YjYl, AiAj = AjAi, YAj = AjY,- L + (6) 
Yr Y r = YfYr, V Y, ),)) . Y, A, A,)) -<),,)). Y, Y, 1 

where Sij = 1 if i = j and otherwise. For convenience we will write (Y i ~) k = Y~ k . 
Let T be the field of functions of n integer arguments a\, a 2 , . . . , a n . The algebra A 
acts on this field 5 , where 

(Yi ■ F)(a 1 , a 2 , • • • , a n ) = F(a ly a;_i, a, + 1, a i+1 , . . . , a n ) , (7) 
(Ai ■ F)(ai, a 2 , • • • , a n ) = aiF(a l7 a 2 , . . . , a n ) . 

Let us turn back to the problem of calculating Feynman integrals. The left-hand 
sides of relations (0) can be represented as elements of the ring A applied to F\ we 

4 An algebra over a field is a vector space over this field and a ring at the same time. 

5 (i) for any a 6 A and / £ T we have an element a ■ f € T\ (it) for any a, b £ A and f,gET 
we have (a + b) ■ (/ + g) = a-f + a- g + b-f + b-g; (Hi) for any a, b € A and / € T we have 
(ab)-f = a-(b-f). 
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will denote these elements by fx, . . . , f n . Now, for F(ax, . . . , a n ) defined by (JTJ), we 
have 

/ 1 .F = 0or(/ i .F)(a 1 ,...,o n ) = (8) 

for all i. Let us generate a (left) ideal X by the elements fx, . . . , f n . We will call X 
the ideal of the IBP relations. Obviously, 

f-F = 0, or {f-F){ ai ,...,a n ) = for any /el. (9) 

Our goal is to express the value of F at an arbitrary point (ax, 0,2, ■ ■ ■ , a n ) in terms 
of the values of F in a few specially chosen points, i.e. master integrals. This problem 
can be solved similarly to the algebraic problem described in Section 2. Consider, for 
example, the case, where all the indices Oj are positive. Then one has 

F( ai , a 2 , . . . , a n ) = (Y^ 1 • • • Y^ 1 ■ F)(l, 1, . . . , 1) . (10) 

The idea of the method is to reduce the monomial Y 1 ai ~ . . . Y® n ~ x modulo the ideal 
of IBP relations. Let us consider a trivial example of such a situation. 
Example 1. One-loop vacuum massive Feynman integrals 

Let us forget that these integrals can be evaluated explicitly, in terms of gamma 
functions. The IBP identity 

A " k W k W^Y^' (12) 

leads to the relation 

(d - 2a + 2)F(a - 1) - 2(a - l)m 2 F(a) = . (13) 

We see that any Feynman integral F(a) where a > 1 can be expressed recursively 
in terms of one integral F(l) = I\ which we therefore qualify as a master integral. 
(Observe that all the integrals with non-positive integer indices are integrals without 
scale and are naturally put to zero within dimensional regularization.) 

Let us demonstrate how the reduction procedure can lead to the same result. 
(We realize that this way is more complicated in this simple situation. However, we 
will see later that its generalization provides simplifications and enables us to solve 
complicated problems.) The IBP relation ()13j) gives us one element / = 2m 2 AY — 
(d — 2 A) G A (the element (/ ■ F)(a - 1) is the left-hand side of flU). Set X = Af 
(for any g G A and F G T we have gf ■ F=0). We have 

2m 2 {A + a- 2)Y a ~ 1 = (2m 2 {A + a - 2)F a ~ 1 - Y a ~ 2 f) + Y a ~ 2 f 

= {2m 2 (A + a - 2)Y a - 1 - Y a - 2 (2m 2 AY - (d - 2A))) + Y a - 2 f 

= Y a - 2 {d - 2A) + Y a - 2 f ={d-2A-2a + 4)Y a ~ 2 + Y a ~ 2 f . (14) 
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The relation 2m 2 (A + a - 2)Y a ~ 1 = (d - 2A - 2a + A)Y a - 2 + X u where X x G X, 
represents one step of the reduction procedure. If we stop the reduction at this point, 

we get 

2m 2 (a - l)F(a) = (2m 2 (a - 1) ■ F)(a) = (2m 2 (A + a - 2)Y a - 1 ■ F)(l) 
= {{d -2A-2a + A)Y a - 2 ) ■ F)(l) + (Y a - 2 f ■ F)(l) 

= (d - 2a + 2)F(a - 1) + ((A - l)Y a - 2 • F)(l) = (d - 2a + 2)F(a - 1) , (15) 

i.e. the equation (|T3|) we started from. But moving further with the reduction modulo 
X we obtain 

P 1 (A,a,m)Y a - 1 = P 2 (A,a,m) + X' , (16) 

where P\ and P 2 are polynomials obtained during the reduction and X' G X (note 
that this algorithm is constructive and is realized as a computer code). Now we can 
apply this equation to F and take the value at 1: 

P 1 (l,a,m)F(a) = (P 1 (A,a,m)Y a - 1 ■ F)(l) 

= ((P 2 (A, a, m) + X') ■ F) (1) = P 2 (l, a, m)F(l) . (17) 

It is enough to notice that Pi is a product of the leading coefficients (the formal 
definition in the non-commutative case will be given later) of /, Yf and so on. Thus 
Pi(l,a, m) is the product of the leading coefficient 2m 2 A of / with A replaced with 
all integers from 1 to a — 1, hence non-zero. After dividing by this value we obtain 
the needed representation. 

Thus it looks tempting to generalize the standard reduction procedure and re- 
duce the monomial Y^~ l . . . Y£ n ~ l so that the resulting polynomial has a smaller 
degree in a certain sense. In this case we would represent F(a\, a 2 , ... , a n ) as a linear 
combination of F(a' v a' 2 , . . . , a' n ) for 'smaller' a[. 

This method works indeed, but first we need to introduce some notation. We will 
say that an element X G A is written in the proper form if it is represented as 

i=^ Ci (4..,4)n^. as) 

i 

where Cj are polynomials and dij are integers. (So, all the operators Ai are placed 
on the left from the operators 5^.) Obviously any element X G A has a unique 
proper form. We will say that an element of A is a monomial if in its proper form 
only one coefficient function Cj is non-zero. We will say that the degree of a monomial 
c(A h . . . , A n ) Hi Y l dl is {di, . . . , d n }. We will say that a monomial c(A u . . . , A n ) Hi Yf* 

is divisible by a monomial c(A%, . . . , A n ) [[ i Y { 1 if d\ > di for all i. 

Let us define a subalgebra A + C A generated as an algebra by Yj and A, (but 
not Y^ 1 ) and set X + = X fl A + . Obviously, X + is an ideal in A + and AX + = 
X. In the same way as in the classical situation, we introduce the notion of an 
ordering, leading term, highest degree of an element of A + and the reduction modulo 
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an ideal. The only problem is that the leading coefficient of an element of A + is now 
a polynomial function, so it does not generally have an inverse element. Thus the 
reduction procedures lead us to a relation 

c (A u . . . , A^Y?- 1 . . . Y^- 1 = c Mii ■ ■ ■ > K) II y f"' + X', (19) 

3 1 

where X' e X, q are polynomials in A,-, and none of the monomials on the right-hand 
side of the relation is divisible by a leading monomial of an element of the basis X + . 
Applying this equality to F and taking the value at (1, . . . , 1) we derive 

q F(a 1 , . . . , a n ) = ^ qjF(d h j + 1, . . . , d nJ + 1) , (20) 

3 

where do not depend on aj. In the case where q is non-zero we can divide by it 
and obtain the desired representation. 

As in the classical situation, one says that a finite set {fi, . . . , /&} is a Grobner 
basis of an ideal X if any element / G X is reduced by this procedure to zero. The 
number of different degrees {di,...,d n } arising on the right-hand side of (|T9"j) is 
minimal possible if we have a Grobner basis. A Grobner basis for the ideal of IBP 
relations can be constructed using the Buchberger algorithm (a generalization of the 
algorithm explained in Section 2). 

Let us now illustrate these points using a very simple 

Example 2. The family of one-loop massless propagator integrals 

FKa2) = / (PM( d g-^ )a2 - (21) 

We have the boundary conditions, F(ai, 02) = if a\ < or a 2 < 0, which correspond 
to putting to zero any integral without scale within dimensional regularization. As it 
is well known, this integral can be evaluated explicitly: 

p/ , , m (~l) ai+a2 r(ai + a 2 + e - 2)T(2 - e - qQrg - e - a 2 ) 

° 2) = (_ g2)ai+ . 2 - 2+£ r(a 1 )r(a 2 )r(4-a 1 -a 2 -2.) ' (22) 

but let us forget about this and consider the problem of the reduction to master 
integrals. 

The two IBP identities 

ddk dk (V 2 ) a i((g-fc) 2 )«* = °) ' (23) 
with I = k and I = q give the following two IBP relations 

d - 2oi - a 2 - a 2 2+(l~ - g 2 ) = , (24) 
a 2 -ai-a 1 l + (g 2 -2")-a 2 2+(l"-g 2 ) = 0. (25) 
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These relations are defined by the elements 

fl = d ^2a l -a 2 -a 2 Y 2 (Y- 1 -q 2 ) , (26) 
h = (a 2 -a 1 )Y 1 Y 2 -a 1 Y 1 (q 2 -Y 2 - 1 )-a 2 Y 2 (Y 1 - 1 -q 2 ). (27) 

which generate the ideal of IBP relations. Let us multiply (|2fijl by Y\ and ()27|) by 
Y{Y 2 to obtain a basis of X + G A + : 

f[ = (d-2a 1 -a 2 -2)Y 1 -a 2 Y 2 (l-q 2 Y 1 ) , (28) 
f 2 = a2 - ai -(a 1 + l)Y 2 (q 2 Y 2 -l) + (a 2 + l)Y 2 2 (q 2 Y 1 -l) . (29) 

If we introduce an ordering for polynomials in the operators and define the 
corresponding reduction procedure modulo the operators f[ and f 2 we shall obtain 
the possibility to represent any given monomial as 

Y^-%^ 1 = rj[ + r 2 f 2 + Y, CijYt^t 1 , (30) 

where r 1 and r 2 are some elements of the ring A. So, if we act by this relation on F 
and take the value at (1,1) we shall obtain 

F(a u a 2 ) = Y°ij F (hj) ■ ( 31 ) 

We discover, however, that the set of integrals that appear on the right-hand side 
of these relations obtained for various a\ and a 2 is infinite. The Buchberger algorithm 
described in Section 2 leads, with the degree-lexicographic order, to the Grobner basis 
consisting of the following two elements: 

gi = 2a 1 Y 1 - da i Y 1 + 2a\Y x - 2a 2 Y 2 + da 2 Y 2 - 2a 2 2 Y 2 (32) 
g 2 = 4aiYi - 4daiYi + d 2 a 1 Y 1 + 8a? Yi - Ada\Yi + Aa\Yi + 2a x a 2 Y 1 
—da\a 2 Y\ + 2a\a 2 Yi + 2a\a 2 Y 2 — da\a 2 Y 2 + 2a 2 a 2 Y 2 — 4q 2 a 2 Y 2 
+dq 2 a 2 Y 2 - Qq 2 a 2 2 Y 2 2 + dq 2 a 2 2 Y 2 - 2q 2 a\Y 2 . (33) 

Now, the reduction modulo these two elements provides only a finite number of 
integrals in the corresponding relations (|31)l . In fact, the degree of g\ is (1, 0) and the 
degree of g 2 is (0, 2), so we meet just the two integrals in this set, F(l, 1) and F(l, 2), 
and call them master integrals. For example, we have 

F(2,3) J < '-^- 5) F(l l 2). (34) 

However, we do not obtain a connection of F(l, 1) and F(l,2), although we know, 
due to explicit solutions of the reduction procedure, that they are connected. This is 
a disturbing point. Of course, it is preferable to have only one master integral in this 
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trivial example, so that we are going to develop an algorithm which reveals a minimal 
number of the master integrals at least in simple examples. 

Starting from example 3, we will have at least one more complication: the variables 
a,i, generally, can be not only positive but also negative. (In the previous example, a\ 
and ci2 were positive due to the boundary conditions.) Generally, we have to consider 
each variable a, to be either positive or non-positive. Of course, for every family of 
Feynman integrals, there will be some boundary conditions. (In particular, if all the 
arguments a, are non-positive any Feynman integral is zero.) 

Thus, if we have a family of Feynman integrals, F(a±, . . . ,a n ), we are going to 
consider each variable Oj to satisfy a« > or < 0. Consequently, we have to consider 
2 n regions that we shall call sectors and label them by subsets v C {1, . . . , n}. The 
corresponding sector a v is defined as {(a 1; . . . , a n ) : > if % G v , aj < if i G" 

In the sector where all are positive, we considered the ring A + C A and the 
operators Yi as basic operators. (See the previous example.) Quite similarly, in a given 
sector <j v it is natural to consider the subalgebra A v C A generated by the operators 
Ai and the operators Yi for i e v and Y~ x for other i. Within this definition we have 

A {l,...,n} = 

Thus the first idea is to construct a Grobner basis for each of the 2 n sectors, or 
at least for all non-trivial sectors. (We call a sector trivial if all the given Feynman 
integrals are identically zero in it due to boundary conditions.) This approach however 
faces many problems: 

1. Each of the non-trivial sectors will give us at least one point where we have to 
evaluate F; 

2. The number of points where we have to evaluate F in a given sector is generally 
greater than the real number of master integrals (In the last example there is 
only one master integral but we obtain F(l, 1) and F(l, 2) after the reduction); 

3. Even if one has constructed all the needed bases, the reduction may fail in cases 
where the coefficient go i n ec L- f)20|) is zero. This problem arises because all 
leading coefficients are polynomials in Ai, so that they can be equal to zero at 
certain points; 

4. Although the method leads us theoretically to constructing a Grobner basis, all 
known practical implementations fail to work even already in four-dimensional 
examples. 

Therefore this specialization of the Buchberger algorithm turns out to be com- 
pletely impractical in sufficiently complicated examples. Our algorithm is a certain 
modification of the Buchberger algorithm. To characterize it we need to introduce 
some notation. 
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Let A^ = ® v iqvA u ' C A. First of all let us define a sector- reduction, or s- 
reduction of an element / G A^ modulo a basis of the ideal X u = X n *4". Take the 
proper form of / and let f v be the sum of the terms in this decomposition that lie in 
A u . If /" is equal to zero the s- reduction stops. Otherwise we look for a monomial 
g G A^ and a coefficient c G /C such that the degrees of / and gfc for some element 
of the bases fi coincide, that (cf — gfi) u is zero or its degree is smaller and that the 
value of c at the point (ai, . . . , a n ) is non-zero, where a, = 1 if i G v and otherwise. 
The procedure is repeated while possible. 

A sector basis, or, an s-basis for a sector a v is a basis of the ideal l v such that 
the number of possible degrees f u , where / is the result of the s-reduction, is finite. 
Such a basis provides the possibility of a reduction to master integrals and integrals 
whose indices lie in lower sectors, i.e. a u > for v 1 C v. 

To prove that an s-basis always exists is an open problem, but in all our examples 
they do exist, and it turns out that in all known examples where one can construct a 
Grobner basis it is an s-basis as well, although this does not follow from the definition. 

If v — then an s-basis is a Grobner basis but generally it is not. Since the 
sector <70 is trivial, we do not have to construct a single Grobner bases. Still, it is 
most complicated to construct s-bases for minimal sectors (a sector a v is said to be 
minimal if it is non-trivial but all lower sectors are trivial). 

Having constructed s-bases for all non-trivial sectors we have an algorithm to 
evaluate F at any point. Indeed, we choose a sector containing the point we need, 
run the s-reduction algorithm for this sector, expressing F in terms of some master 
integrals and values for lowers sectors, then repeat the procedure for all those sectors. 
Eventually we reduce F to the master integrals. 

The Buchberger algorithm leads us to constructing a Grobner basis that is hope- 
fully an s-basis, but this has no use for us since this does not simplify anything. The 
second important point is that the Buchberger algorithm can be terminated when the 
Grobner basis is not yet constructed but the 'current' basis already provides us the 
s-reduction, so that it is an s-basis (we have criteria that show whether a basis is an 
s-basis) . 

Let us illustrate how this idea works on the same example. The initial basis turns 
out to be an s-basis. First, let us observe that the degree of f[ is (1,1) and the 
leading coefficient is q 2 a2, i.e. is a non-zero function in the positive sector, hence 
we are capable of making reduction steps if the highest degree of an element being 
reduced is different from (/,0) and (0,/). Now let / be a polynomial whose highest 
degree is (1,0) and the leading coefficient is c. Then 

/' = q \l - i + ai )f + cYt 2 Y 2 - l f 2 (35) 

is an element of A. Let us take the proper form of / (implying that the numbers d it j 
can be now negative) and calculate its highest degree without paying attention to the 
terms with negative dij. Obviously it is smaller than the degree of /. Now if we take 
the value of /' at (1, 1) we will have the elements like F(j, 0) among others. But the 
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boundary conditions state that they are equal to zero, so that we have nothing to 
worry about. Now we can move further and reduce /' as we did before. Finally we 
get 

f = X + c(A 1 ,A 2 )+Y (36) 

where X e X, c is a rational function and Y is an element of A such that in its proper 
form there are no degrees di,d 2 with d\ > and d 2 > 0. Taking the value at (1, 1) 
we obtain the reduction of any value of F to the value of the master integral F(l, 1). 
(See also [13].) 



4 Examples 

Let us consider a modification of Example 2; now we have a non-zero mass mi = m: 
Example 3. Propagator integrals with the masses m and 0, 

/d d k 
(jfc2_ m 2)a l[(g _ fc) 2 ]fla " ( 3? ) 

The integrals are zero if a% < 0. The corresponding IBP relations generate the 
following elements: 

/i = d — 2a\ — a 2 — 2m 2 a\Y\ — m 2 a 2 Y 2 + q 2 a 2 Y 2 — a 2 Y 2 Y-^ 1 

f 2 = a 2 — a\ — m 2 aiYi - q 2 a 1 Y 1 - m 2 a 2 Y 2 + q 2 a 2 Y 2 - a 2 Y 2 Y{ x + a{YiY 2 1 . 

We have to consider two sectors, <?{i, 2 } and a^. 

Using the lexicographical ordering, we obtain, for the sector cr{i l2 \, the s-basis 
consisting of two elements: 

gn = Y 2 + ai Y 2 + W X Y 2 - dY{Y 2 + a x Y x Y 2 + 2a 2 Y 1 Y 2 + m 2 Y 2 Y 2 
-q 2 Y 2 Y 2 + m 2 a 1 Y 1 2 Y 2 - q 2 a l Y 2 Y 2 , 



gi2 = -3Y X Y 2 + dY{Y 2 - 2a 1 Y 1 Y 2 - a 2 Y x Y 2 - 2m 2 Y 2 Y 2 - 2m 2 a x Y 2 Y 2 - Y 2 



-a 2 Y 2 - m^Y 2 + q 2 Y x Y 2 - m 2 a 2 Y 1 Y 2 + q 2 a 2 Y 1 Y 2 . 
For the sector cr{i}, we obtain the following s-basis: 

g 21 = 1 - a 2 + m 2 Y x - q 2 Y x - m 2 a 2 Y 1 + q 2 a 2 Y 1 - Y 1 Y 2 ' 1 + dY x Y 2 x - 2a 1 Y 1 Y 2 ~ 1 

-a 2 Y x Y 2 x - 2m 2 Y 2 Y 2 1 - 2m 2 aiY 2 Y 2 l , 
g 22 = -2m 2 + 2m 2 a 2 - 2m 4 Fi + 2m 2 q 2 Y 1 + 2m 4 a 2 Y 1 - 2m 2 q 2 a 2 Y 1 - 2Y 2 X + a 2 Y 2 x 

+2m 2 Y 1 Y 2 1 + 2q 2 Y 1 Y 2 1 + 2m 2 a 1 Y l Y 2 1 - m 2 a 2 Y 1 Y 2 1 - q 2 a 2 Y 1 Y 2 1 

+2m A Y 2 Y 2 1 + 2m 2 q 2 Y 2 Y 2 1 + 2m i a l Y 2 Y 2 1 + 2m 2 q 2 a l Y 2 Y 2 1 

—dY{Y 2 2 + 2a l Y l Y 2 2 + a 2 Y x Y 2 2 . 
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Figure 1: Two-loop propagator diagram 



The reduction based on the two constructed s-sectors reveals two master integrals, 
F(l, 1) and F(1,0), in accordance with results obtained by other ways. (See, e.g., 
Chapters 5 and 6 of [T5] and [T5].) 

Example 4- Two-loop massless propagator diagram of Fig. 1. 

The corresponding family of Feynman integrals is 

f f d d kd d l 
F(a u a 2 , a 3 , a 4) a 5 ) - J J (F)fll[(g _ fc)2]fl2(/2)a3[(g _ /)2]a4[(fc _ /)2]a5 • (38) 

There are boundary conditions which correspond to setting to zero integrals without 
scale: F(a\, a 2 , a 3 , a 4 , a 5 ) = , if aj, as < for i — 1, . . . , 4, or a 1; a 2 < 0, or a 3 , a 4 < 0, 
or ax, a 3 < 0, or a 2 , a 4 < 0. The integrals are symmetrical: 

F(ai, a 2 , a 3 , a 4 , a 5 ) = F(a 2 , a 4 , a 4 , a 3 , a 5 ) = -F(a 3 , a 4 , a 1; a 2 , a 5 ) . 

The corresponding IBP relations generate the following elements: 



fx 


= (d- 


2a x - 


- a 2 


- a 5 ) 


+ a 2 Y 2 (q 2 - 


^f 1 ) 


-a.Y^Y,- 1 - 




h 


= (d- 


a 2 - 


2a 3 


- as) 


+ a 4 F 4 (g 2 - 


IT 1 ) 


- a.Y.iYi 1 - 




h 


= (d- 


d\ — 


a 2 - 


-2a 5 ) 


+ a 1 Y 1 (Y 3 ~ 1 




" 1 ) + a 2 F 2 (y 4 - 1 




h 


= (d- 


a 3 - 


a 4 - 


-2a 5 ) 


+ a 3 Y 3 (Y 1 - 1 




- 1 ) + a A Y A (Y 2 l 




h 


= (d- 


ai — 


2a 2 


- «5) 


+ a 1 Y 1 (q 2 - 




- a 5 y 5 (y 2 - x - 




h 


= (d- 


a 3 - 


2a 4 


- a 5 ) 


+ a 3 Y 3 {q 2 - 


vr 1 ) 


- a 5 y 5 (yr 1 - 


n- 1 ) • 



As is well known, any integral of this class can be reduced, due to IBP relations, in 
a very simple way, to integrals where at least one of the indices is non-positive. Such 
integrals can be evaluated recursively in terms of gamma functions using the one-loop 
integration formula (|22|). Let us point out that physicists often stop the reduction 
whenever they arrive at integrals expressed in terms of gamma functions. Imagine, 
however, that we want to know the whole solution of the reduction procedure, i.e. a 
reduction to a minimal number of the master integrals. Then this example turns out 
be not so trivial and provides a good possibility to test our algorithms. 
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This can be done in various ways. In our approach, we apply our algorithm to 
construct s-bases corresponding to the sectors 0"{i,2,3,4,5}) °"{2,3,4,5} as well as three 
more (symmetrical) sectors, the minimal sector <7{i i2;3i 4}, the minimal sector 0"{2, 3 ,5} 
as well as one more (symmetrical) sector. In the first three cases, we used the degree- 
lexicographical ordering, and in the last case, some special ordering. 

For example, let us present the s-basis associated with 0"{i 23 4}: 

91 = 2Y 1 Y 2 Y 3 Y 5 ~ 1 - a 5 Y 1 Y 2 Y 3 Y 5 - 1 - 2Y 1 Y 2 Y 4 Y 5 ~ 1 + a 5 YiY 2 Y 4 Yf 1 - 2Y 1 Y 3 Y 4 Yf 1 

+a 5 Y 1 Y 3 Y 4 Y 5 ' 1 + 2Y 2 Y 3 Y 4 Y 5 ' 1 - a h Y 2 Y 3 Y 4 Y^ x + 2q 2 Y x Y 2 Y 3 Y 4 Y^ x 
—q 2 o,5YiY 2 Y 3 Y 4 Yf 1 + 2q 2 Y l Y 2 Y 2 Y^ 1 - q 2 a 5 Y 1 Y 2 Y 4 Y 5 ~ l - 2q 2 Y 1 Y 3 Y 2 Yf 1 
+q 2 a 5 Y 1 Y 3 Y 2 Y 5 - 1 - 2q 2 Y 2 Y 3 Y 2 Y,- 1 + q 2 a 5 Y 2 Y 3 Y 2 Y^ - (q^Y^Y^Yf 1 
+d{q 2 ) 2 Y l Y 2 Y 3 Y 2 Y,- 1 - (q 2 ) 2 ■a 3 Y 1 Y 2 Y 3 Y 4 2 Yf 1 - {q 2 ) 2 a 4 Y x Y 2 Y 3 Y 2 4 Y^ 

-2( g 2 ) 2 a 5 y 1 y 2 y 3 n 2 n _1 + {q 2 ) 2 Y 2 YiY 2 Y^ + {q 2 ) 2 a,Y 2 Y 2 Y 2 Y^ 

+2(q 2 ) 2 Y 1 Y 3 Y?Y 5 - 1 + (q 2 f 'a^Y^Yf 1 - Y,Y 2 Y 2 Y^ 2 - a^Y^Yf 2 

+a 3 yiy 2 y3y4y5~ 2 - a 4 Y x Y 2 Y 3 Y 4 Y^ 2 + yy^yf 2 + a 4 Y x Y 2 YiY^ 2 
-6g 2 y 1 y 2 y 3 y 4 2 y 5 - 2 + 2rfg 2 y 1 y 2 y 3 y 4 2 y 5 - 2 - 3g 2 a 3 y 1 y 2 y 3 y 4 2 y 5 - 2 
^Vyyyyfy,- 2 - 2 g 2 a 5 yiy2y 3 y 4 2 yr 2 - 2g 2 y 1 y 2 y 4 3 y 5 - 2 
-gVyyyfyr 2 , 

92 = Y{Y 2 Y 2 - a 5 YiY 2 Y 2 - Y 2 Y 3 Y 2 + a 5 Y 2 Y 3 Y 2 - SY^Y^Yf 1 

+dY 1 Y 2 Y 3 Y 2 Y 5 - 1 - 2a 3 Y 1 Y 2 Y 3 Y 2 Y 5 - 1 - a 4 Y x Y 2 Y 3 Y 2 4 Y^ x - a^Y^Y^Yf 1 
-2Y x Y 2 YlY^ x - a 4 Y x Y 2 YlY^ x + 2g 2 y 1 y 2 y 3 y 4 3 y ) - 1 + q 2 a 4 Y 1 Y 2 Y 3 Y 4 3 Y 5 ~ 1 , 

93 = -yyyy + ^y x y 2 y 3 y 4 + yy 2 y 4 2 - asyyy 2 + y x y 3 y 2 - ^y x y 3 yi 

-Y 2 Y 3 Y 2 + a 5 Y 2 Y 3 Y 2 + Y^Y^Y' 1 + a^Y^Y.Y,' 1 + Y^Y^Y' 1 
-a 3 Y x Y 2 Y 3 YlY^ x + a 4 Y x Y 2 Y 3 YlY^ - q^Y^Y^ 1 - q 2 a^Y.Y^Yf 1 
-2Y 1 Y 2 Y 4 %- 1 - a 4 Y x Y 2 Y 3 4 Y^ + 2g 2 y 1 y 2 y 3 y 4 3 y- 1 + q 2 a 4 Y 1 Y 2 Y 3 Y 4 3 Y 5 ~ 1 , 
g 4 = -Y,Y 2 Y 2 + a b Y,Y 2 Y 2 + Y 2 Y 3 Y 2 - a 5 Y 2 Y 3 Y 2 - 2Y 1 Y 2 Y 3 Y 4 %- 1 

+dY 1 Y 2 Y 3 Y 2 Yf 1 - 2a 1 Y 1 Y 2 Y 3 Y 2 Y 5 - 1 - a 2 Y x Y 2 Y 3 Y 2 4 Y^ - a^Y^Y^ 1 

-Y 2 Y 3 Y 2 Y b ~ l - a 2 Y 2 2 Y 3 YiY^ x + q 2 Y x Y 2 Y 3 Y 2 4 Y^ 1 + q 2 a 2 Y x Y 2 Y 3 Y 2 4 Y.; x , 

g 5 = Y ± Y 2 Y 3 Y 4 - a 5 Y,Y 2 Y 3 Y 4 - Y 1 Y 2 Y? + a 5 Y,Y 2 Y 2 - Y 1 Y 3 Y 2 + a^Y^Y 2 

+Y 2 Y 3 Y 2 - a 5 Y 2 Y 3 Y 2 + YfY^Y^ 1 + a^Y^Yf 1 - a.Y^Y^Yf 1 
+a 2 Y 1 Y 2 Y 3 Y 2 Y 5 - 1 - q 2 Y 2 Y 2 Y 3 Y 2 Y^ 1 - q 2 ' ai Y 2 Y 2 Y 3 Y 2 Y^ - Y 2 2 Y 3 Y 2 Y 5 ~ l 
-a 2 Y 2 Y 3 YlY^ x + q 2 Y x Y 2 2 Y 3 YlY^ x + q 2 a 2 Y x Y 2 2 Y 3 Y 2 Y^ , 

g 6 = -q 2 Y x Y 2 Y 2 + q 2 a 5 Y 1 Y 2 Y 2 + q 2 Y 2 Y 3 Y 2 - q 2 a 5 Y 2 Y 3 Y 2 + 2Y l Y 2 Y 3 Y 4 Y^ 1 

-a^Y x Y 2 Y 3 Y 4 Y^ - 2Y X Y 2 Y 2 Y^ + a^Y^Yf 1 - 2Y 1 Y 3 Y 2 Y 5 ~ 1 + a^Y^Y^ 1 
+2Y 2 Y 3 Y 2 Y 5 ~ 1 - a b Y 2 Y 3 Y 2 Y^ + 2g 2 y 1 y 2 y 3 y 4 2 y ) - 1 + q 2 a 3 Y 1 Y 2 Y 3 Y 2 Y 5 - 1 
-q 2 a h Y x Y 2 Y 3 YlY^ + q 2 Y 2 Y 2 Y 2 Y^ + q 2 a 3 Y 2 Y 2 Y 4 %- 1 + 2q 2 Y 1 Y 2 Y 4 %- 1 
+q 2 a 4 Y 1 Y 2 Y 4 %- 1 + 2g 2 y 1 y 3 y 4 3 yr 1 + q 2 a^Y^Yf 1 - 2(q 2 ) 2 Y 1 Y 2 Y 3 Y 4 %- 1 
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-{q 2 ) 2 a 4 Y 1 Y 2 Y 3 YiYf 1 - Y 1 Y 2 Y 2 Y 4 Y<f 2 - a 3 Y 1 Y 2 Y 3 2 Y 4 Y 5 - 2 - Y^Y^Yf 2 
+a 3 Y 1 Y 2 Y 3 Y 2 Y 5 - 2 - a±Y x Y 2 Y 3 Y 2 Y^ + 2Y 1 Y 2 Y 4 3 Yf 2 + a^Y^Y' 2 
—4q 2 YiY 2 Y 3 Y 4 3 Yf 2 - 2q 2 a 4 Y 1 Y 2 Y 3 Y 3 Y 5 ~ 2 , 
g 7 = q 2 Y,Y 2 Y 2 - q 2 a,Y l Y 2 Y 2 - q 2 Y 2 Y 3 Y 2 + q 2 a b Y 2 Y 3 Y 2 - 2Y 1 Y 2 Y 3 Y 4 Y 5 ~ 1 

+a 5 Y 1 Y 2 Y 3 Y i Y 5 - 1 + q 2 Y x Y 2 2 Y 3 Y 4 Y^ + q 2 a 2 Y 1 Y 2 Y- 3 Y^ x + 2Y 1 Y 2 Y 2 Yf 1 
-a 5 Y x Y 2 Y 2 Y^ + q 2 Y 2 Y 2 Y 2 Y- 1 + q 2 aiY 2 Y 2 Y 2 Y" 5 ~ 1 + 2Y 1 Y 3 Y 2 Y ; r 1 
-a^YzYlYf 1 - 2Y 2 Y 3 Y 2 Yi 1 + a^Y 2 Y 3 Y 2 Y^ x + 2q 2 Y l Y 2 Y 3 Y 2 Y^ 1 
+q 2 a 1 Y 1 Y 2 Y 3 Y 2 Y 5 ~ 1 - q 2 a^YiY 2 Y 3 Y 2 Y^ 1 + q 2 Y 2 Y 3 Y 2 Y^ + q 2 a 2 Y 2 Y 3 Y 2 Y^ 
~{q 2 ) 2 Y x Y 2 2 Y 3 YlY^ - {q 2 ) 2 a 2 Y x Y 2 Y 3 YlY^ - Y 2 Y 3 Y 2 Y 5 ~ 2 - ai Y 2 Y 3 Y 4 %~ 2 
+a 1 Y 1 Y 2 Y 3 Y 2 Yf 2 - a 2 YiY 2 Y 3 Y 2 Y?~ 2 + Y 2 Y 3 Y 2 Y<r 2 
+a 2 Y 2 Y 3 Y 2 Y,- 2 - 2q 2 Y 1 Y 2 2 Y 3 Y 2 Y 5 ~ 2 - 2q 2 a 2 Y 1 Y 2 2 Y 3 Y 2 Y,- 2 . 

The reduction based on the constructed s-sectors reveals three master integrals, 
F(l, 1, 1, 1, 0), F(0, 1, 1, 0, 1) and F(l, 0, 0, 1, 1) (the last two of them are equal be- 
cause of the symmetry), in accordance with results obtained by other ways. (See, 
e.g., Chapters 5 and 6 of [T5].) 
Our last example is 

Example 5. Two-loop Feynman integrals for the heavy quark static potential 
corresponding to Fig. 2 with v ■ q = 0. 




Figure 2: Feynman diagram contributing to the three-loop static quark potential. A 
wavy line denotes a propagator for the static source and the dotted line denotes the 
scalar propagator with the index shifted by e. 

We shall consider diagrams contributing to the three-loop static potential corre- 
sponding to Fig. 2. They are obtained from the corresponding two-loop diagrams by 
inserting a one-loop diagram into the central line. Indeed, the integration over the 
loop-momentum of the insertion can be performed explicitly, by means of (|22|). and 
one obtains, up to a factor expressed in terms of gamma functions, Feynman integrals 
of Fig. 2, where the index of the central line 6 is + e = as + (4 — d)/2 with integer 

6 A more general case, with 0,5 — > 05 + re and integer r, relevant to r-loop massless insertions can 
be considered on the same footing. 
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a§. So, we arrive at the following family of integrals: 7 



F(ai, a 2 , a 3 , a 4 , a 5 , a 6 , a 7 ) — 




d d fcd d Z 



(_Jfe2)ai[_(jfe _ g )2]a 2 (_/2)a 3 [_(/ _ g )2] 



1 



(39) 



x 



[-(k - /) 2 ] a 5+ £ (-w • k) a *(-v ■ Z) a7 ' 



We have turned to the —k 2 dependence of the propagators because this choice is more 
natural when at least one index, as + e is not integer. The integrals are symmetrical: 



They are equal to zero, if ax, a 3 < 0, or 02, a 4 < 0, or ax, a 2 , a 6 < 0, or a 3 , a 4 , a 7 < 0. 

The IBP relations generate the following elements in the case where a 5 is not 
shifted by e (i.e. the case of the diagrams relevant to the two-loop static quark 
potential considered in [Til IS]) 



h = ( d _2ax-a 2 -a 5 -a 6 )-a 2 Y 2 (q 2 + Y 1 - 1 )-a 5 Y 5 (Y 1 - 1 -Y 3 - 1 ), 

h = (d-a 2 -2a 3 -a 5 -a 7 )-a 4 Y 4 (q 2 + Y 3 - 1 )-a 5 Y 5 (Yi 1 -Y 1 - 1 ), 

h = (d-a 1 -a 2 -2a 5 -a e )+a 1 Yx(Y 3 - 1 -Y 5 - 1 ) + a 2 Y 2 (Y 4 - 1 ~Yf 1 ) + a 6 Y 6 Y 7 - 1 , 

h = (d-a 3 -a 4 -2a 5 -a 7 )+a 3 Y 3 (Y 1 - 1 -Yi 1 ) + a 4 Y 4 (Y 2 - 1 -Y 5 - 1 ) + a 7 Y 6 - 1 Y 7 , 

f 5 = (d-a 1 -2a 2 -a 5 -a 6 )-a 1 Y 1 (q 2 + Y 2 ~ 1 )-a 5 Y 5 (Y 2 - 1 -Y 4 ~ 1 ), 

h = (d-a 3 -2a 4 -a 5 -a 7 )-a 3 Y 3 (q 2 + Y 4 ~ 1 )-a 5 Y 5 (Y 4 - 1 -Y 2 - 1 ), 

ft = 2a 1 Y 1 Y 6 - 1 + 2a 2 Y 2 Y 6 ~ 1 + a 5 Y 5 (Y G - 1 -Y 7 - 1 )-v 2 a 6 Y 6 , 

f 8 = 2a 3 Y 3 Y 7 - 1 + 2a 4 Y 4 Y 7 - 1 -a 5 Y 5 (Y 6 - 1 -Y 7 - 1 )-v 2 a 7 Y 7 . 



So, the IBP elements we need are obtained from these by replacing 05 with 05 + e. 

Our algorithm works successfully in this example and gives us a family of s-bases 
which provide the possibility of a reduction to master integrals. The elements of 
the bases are rather lengthy, typically, with hundreds of terms, so that we do not 
present them in this short paper. These s-bases correspond to the following sec- 
tors: &{x,2,3A'3fi,7}-> <T {2,3,4,5,6,7}, 2,3,4,5, 7} > °"{3,4,5,6,7} j °"{2,3,5,6,7}> °"{2,3,4,5,7} > °"{2,3,4,5,6}> 

2,3,4,5}, o"{2,3,4,5}, ^"{2,3,5,6} and other sectors obtained by the symmetry transforma- 
tions. 



We obtain the master integrals: h = F(l, 1, 1, 1, 0, 1, 1), I 2l = F(l, 1, 1, 1, 0, 0, 1), 
I 22 = F(l, 1, 1, 1, 0, 1, 0), I 3 = F(l, 1,1,1, 0, 0, 0) . We have I 2X = I 12 = I 2 because of 
the symmetry. We also obtain J 51 = F(l, 0, 0, 1, 1, 1, 1), 7 71 = F(l, 0, 0, 1, 1, 0, 1), I 81 = 
F(1,0,0, 1,1,1,0), I 4l = F(l, 0,0, 1,1,0,0). We have J 71 = I 81 = h because of the 
symmetry. Moreover, we have other copies, 152,^72,-^82,-^42, of this last family of the 
master integrals which are obtained by the symmetry transformation (1 <-> 2,3 «-> 



7 These integrals with the integer index of the central line contribute to the two-loop static quark 
potential and were calculated in |14ll8). 



F(ax, a 2 , a 3 , a 4 , 05, a 6 , a 7 ) — F(a 2 , ax, a 4 , a 3 , as, a§, a 7 ) — F(a 3 , a 4 , a 4 , a 2 , as, a 7 , a§) . 
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4). We also obtain J 61 = F(0, 0, 1, 1, 1, 1, 0), I el = F(0, 0, 1, 1, 1, 2, 0) as well as the 
corresponding symmetrical family. 

To calculate the master integrals one can use the threefold Mellin-Barnes repre- 
sentation 



F(di, <2 2 , «3, Gt4, 05, Qq, a-j) 



2*r-i(v 2 ) 



-0167/2 



X 



X 



(2ttz) 



+JOO 



+JOO 



dzidz 2 dz; 



1^=3,4,5,7 r(a/)r(4 - a 3 457 - 2e)(Q 2 )<™-4 + 2 £+ a 67 /2 

r(a 12345 + a 67 /2 + 2e-4 + ,2 3 ) 



r(a 345 + a 67 /2 + £ - 3/2 + 21 + z 2 + z 3 ) 



r(a 3 + zi + -2 3 )r(a 4 + z 2 + ^ 3 )r(a 345 + a 7 /2 + s - 2 + zi + z 2 + z 3 ) 



r(oi - zi)r(aa - ^ 2 )r(8 - ai 267 - 2a 345 - 4e - z x - z 2 - 2z 3 ) 
xr(a 345 + a 7 /2 + e - 3/2 + z x + z 2 + z 3 )r(4 - 2a 34 - a 57 - 2e - z y - z 2 - 2z 3 ) 
xT(4 - ai 345 - a 67 /2 - 2e - z 2 - z 3 )T(2 - a 345 - e - z x - z 2 - z 3 ) 
xr(4 - a 2345 - a 67 /2 -2e-zy- z 3 )V(-Zy)T(-z 2 )V(-z 3 ) , (40) 

where a 4234 5 = a\ + a 2 + a 3 + a 4 + 05 etc. The integrations over the variables Z{ go 
from —00 to +00 in the complex plane. The contours are chosen in the standard way: 
the poles of gamma functions with the — Zi dependence are to the right of the contour 
and the poles with the +Zi dependence are to the left of it. This representation 
can be derived by applying Feynman parameterization to the sub loop integral over /, 
introducing then three MB integrations and, finally, integrating over k. (See Chapter 
4 of ^H] for details of this method, with multiple examples.) 
We obtain the following results for the master integrals: 



8tt2 
~9e 



16tt 2 40C(3) 
"9" + 3 



O(e) 



171 



d/2. 



-7E<n 



Q 



3+6e, 



( 



Q 



2+6e 



7T 



(l7T 



d/2 



)Q 



2-6e 



6C(3) + ( - + 12C(3)J e 
r(l-2e)r(l -s) 2 T(3s- 



-0(s 2 ) 

1) 



d/2 --y E e) 



Q6e~2 v 2 



r(3-4e)r(l + e 
8C(3) 



{i* d/2 ) 2 Q 



2-6e 



4^ 32^_^, 
9e 9 3 v ; 

4 1 - 2e v /^r (3/2 - 3ef r (1 - 2e) T (3e - 1/2) T (4e - 1) 



{i7i d ' 2 ) 2 Q 2 - 6£ - 



l-2e 



T (3 - 6e) T (2e) T (1 + e) v 
v^Fr (1 - 3s) 2 r (1 - 2e) r (3e) T (4e) 
T (2 - 6e) r (1 + e) f (1/2 + 2e) v 2 



(m d / 2 ) 



i2 v/?r (3/2 - 3e) T (1 - 2e) V (1/2 - e) Y (1 - e) V (3e - 1/2) 

g^-^r (2 - 4e) r (2 - 3e) r (1 + 7) 
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h = It, 



where Q = \/—q 2 and v = \f\P'. 

Observe that some of the integrals are expressed explicitly in terms of gamma 
functions for general values of e while results for some other integrals are presented in 
expansion in a Laurent series in e. The depth of this expansion can be made greater 
whenever necessary. 

For example, we obtain the following reductions to master integrals by our algo- 
rithm: 



1,1, 1,1, i,-i) = -,^ v J,-3/,- 8(< '- 3) ^- 7)(1 ^ 46) ^ 

v ' ' ' ' ' ' ; (3d -10) (d-4) 2 (3d- 14)Q 4 



4(3d-ll)(7d-30), 2 

(d-4)(3d-14)(3d-10)Q 2 6 ' K ' 

3d - 14 r 
F(2, 1,1, 1,1, 1,1) = 

4(d - 3) (d - 2) (2d - 7) [3d - 10) (9d - 40) 

~ (d - 5)(d - 4) (2d - ll)(3d - 16)(3d - 14)Q% 2 4 

3(d-4)(4d- 17)(4d- 15) 16(3d - 13)(3d - 11) - 

"-'s - 77T3 7TT7773 T7\7T3 77\7v;-'6 , (42) 



(2(d-5)(2d- 11)Q 6 ° (2d- ll)(3d- 16)(3d- 14)Q 6 ' 

which can be checked straightforwardly, by evaluating these integrals, in expansion 
in e, using the MB representation (|40p. 



5 Conclusion 

We have developed an algorithm which is a generalization of the Buchberger algorithm 
to the reduction problem for Feynman integrals and modified it in such a way that 
it works at the level of modern calculations. We have described the main features of 
the algorithm. For the examples considered, it works rather fast — these are seconds 
of CPU time for Example 4 and minutes for Example 5, both for constructing s-bases 
and reduction to master integrals. In fact, it has turned out that our algorithm works 
successfully even at a higher level, in a reduction problem with nine indices THI. Still 
to perform more sophisticated calculations, further modifications and optimizations 
are needed. One of possible ways to improve the algorithm is to combine its basic 
points with that of algorithms based on Janet bases [T2] . We hope to report on 
our progress in future publications. We also postpone to solve various mathematical 
problems connected with our algorithm. 

Acknowledgments. We are grateful to I.V. Arzhantsev, V.P. Gerdt, Y. Schroder 
and M. Steinhauser for helpful discussions and careful reading of draft versions of 
the paper. The work of A.S. was supported by the Russian Foundation for Basic 
Research through grant 05-01-00988 and by CRDF through grant RM1-2543. The 



18 



work of V.S. was supported by the Russian Foundation for Basic Research through 
grant 05-02-17645. 

References 

[1] G. 't Hooft and M. Veltman, Regularization and renormalization of gauge fields, 
Nucl. Phys. B 44 (1972) 189; C.G. Bollini and J.J. Giambiagi, Dimensional 
renormalization: the number of dimensions as a regularizing Parameter, Nuovo 
Cim. 12 B (1972) 20. 

[2] K.G. Chetyrkin and F.V. Tkachov, Integration by parts: the algorithm to calcu- 
late beta functions in 4 loops, Nucl. Phys. B 192 (1981) 159. 

[3] S. Laporta, High-precision calculation of multi-loop Feynman integrals by differ- 
ence equations, Int. J. Mod. Phys. A 15 (2000) 5087 |hep-ph/0102033| . 

[4] T. Gehrmann and E. Remiddi, Two-loop master integrals for 7* — > 3 jets: the 
planar topologies, Nucl. Phys. B 601 (2001) 248 | hep-ph/0008287| ; T. Gehrmann 
and E. Remiddi, Two-loop master integrals for 7* — > 3 jets: the non-planar 
topologies, Nucl. Phys. B 601 (2001) 287 hep-ph/ 01 011241 . 

[5] C. Anastasiou and A. Lazopoulos, Automatic integral reduction for higher order 
perturbative calculations, JEEP 0407 (2004) 046 |hep-ph/0404258| . 

[6] O.V. Tarasov, A new approach to the momentum expansion of multiloop Feynman 
diagrams Nucl. Phys. B 480 (1996) 397 hep-ph/9606238|; Connection between 
Feynman integrals having different values of the space-time dimension, Phys. 
Rev. D 54 (1996) 6479 |hep-th/ 96060181 . 

[7] P. A. Baikov, Explicit solutions of the 3-loop vacuum integral recurrence relations, 
Phys. Lett. B 385 (1996) 404 |hep-ph /9603267|; Explicit solutions of the multi- 
loop integral recurrence relations and its application, Nucl. Instrum. Methods A 
389 (1997) 347 |[hep-ph/9611449| . 

[8] V.A. Smirnov and M. Steinhauser, Solving recurrence relations for multi-loop 
Feynman integrals, Nucl. Phys. B 672 (2003) 199 | hep-ph/0307088| . 

[9] B. Buchberger and F. Winkler (eds.) Grobner Bases and Applications, (Cam- 
bridge University Press, 1998). 

[10] O.V. Tarasov, Reduction of Feynman graph amplitudes to a minimal set of basic 
integrals, Acta Phys. Polon. B 29 (1998) 2655 | hep-ph/9812250| ; Computation 
of Grobner bases for two-loop propagator type integrals, Nucl. Instrum. Meth. A 
534 (2004) 293 |hep-ph/0403253"l. 



19 



[11] V.P. Gerdt, Grobner bases in perturbative calculations, Nucl. Phys. B (Proc. 
Suppl), 135 (2004) 2320 hep-ph/0501053| . 

[12] V.P. Gerdt, On computation of Grobner bases for linear difference systems 
|math-ph/0509050l . 

[13] V.P. Gerdt and D. Robertz, A Maple package for computing Grobner bases for 
linear recurrence relations cs.SC/0509070 . 

[14] M. Peter, The static quark- antiquark potential in QCD to three loops, Phys. Rev. 
Lett. 78 (1997) 602 |hep-ph/9610209 ; The static potential in QCD: a full two- 
loop calculation, Nucl. Phys. B 501 (1997) 471 |hep-ph/9 70224 5|; Y. Schroder , 
The two-loop static potential, Phys. Lett. B 447 (1999) 321 |hep-ph/9909520| ; 
B.A. Kniehl, A. A. Penin, V.A. Smirnov, and M. Steinhauser, Nonabelian 
a(s) 3 1 (m(q)r 2 ) heavy- quark- antiquark potential, Phys. Rev. D 65 (2002) 091503 
|hep-ph/0106135| . 

[15] V.A. Smirnov, Evaluating Feynman integrals, (STMP 211; Springer, Berlin, Hei- 
delberg, 2004). 

[16] A.G. Grozin, A.V. Smirnov and V.A. Smirnov, to be published. 



20 



